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I N THE AREAS of animal orientation and 
navigation/ directions are measured in 
various ways. Animals are either kept in 
cages and their movements or preferences 
with respect to a reference direction recorded, 
or they are released and at certain points of 
their path the bearings are measured. The 
purpose of a statistical analysis is to establish 
preferred directions and to compare them 
with geographical lines such as directions to 
breeding places, homeward directions, val- 
leys, mountain ranges, shore lines, and die 
like. It is also important to find significant 
differences between experimental and control 
groups. 

A comprehensive account of statistical 
methods in the analysis of directions for biolo- 
gists was published in reference 1. In the 
past five years, there was considerable inter- 
est in this area by both theoretical and ap- 
plied research workers. A host of new results 
and methods is available today. 

This paper explains and illustrates those 
methods that are immediately useful for 
biologists. No attempt is made to review re- 
sults that are of purely theoretical interest no 
matter how important they may be for future 
statistical research. We will restrict ourselves 
to the two-dimensional or circular case and 
omit the three-dimensional or spherical case. 


An exhaustive bibliography is found at 
the end of this paper for the reader whose 
interest goes beyond the rather limited scope 
of this account. Only those publications are 
listed that are not already quoted in Bat- 
schelet, 1965 ( ref. 1 ) . Recent papers of 
basically mathematical interest can be found 
under such names as Ajne, Beran, Downs, 
Maag, Mardia, Rao, Schach, Stephens, and 
Watson. For the spherical case we refer the 
reader to Bingham (1964), Downs (1966, 
1970) , Downs and Gould (1967), Savary 
(1965), Selby (1964), Stephens (1967, 
1969), Watson (1965, 1966, 1967, 1970). 

NOTATIONS 

We will use the same notations as in ref- 
erence 1. Let O'!, «25 • • • 5 be a sample 
of independent angles taken from a certain 
theoretical distribution. Without loss of gen- 
erality we choose the positive x axis as zero 
direction and associate the positive sign with 
a counter-clockwise rotation (see fig. 1 ) . 

It is convenient to introduce unit vec- 
tors, that is, vectors of unit length point- 
ing in the directions given by the angles oji 
(izz- lj-T^n). The tips of these unit vec- 
tors are located on the circumference of the 
unit circle. 


61 




62 


ANIMAL ORIENTATION AND NAVIGATION 



FIGURE 1. drcular sample represented by unit 
vectors. 

In a rectangular ;Tjy-coordinate system 
the components of the unit vectors are 

Xi == cos ctiy yi sin at (t = 1, • • • , n) (1) 

When a preferred direction or a mean 
direction is to be defined, we associate with 
each direction a unit mass which we locate 
at the tip of the corresponding unit vector 
(fig. 2) . Then we calculate the coordinates 
of the center of gravity by 

^ 2co8a„y=--2w = 
i 2 sin a, (2) 

ft 

Denote the center of gravity by C* The 
vector which points to C is called the mean 
vector, whose components are x and y. 

The length of the mean vector is denoted 
by r and may be calculated from 

r=^(3c2 4.f)i/2 (3) 

An alternative way of calculating r uses 
vector addition. Consider again the unit vec- 
tors with components cos at, sin at (i = 1, 
• • • , n) . The sum of these vectors, also called 
the resultant vector, has components 



FIGURE 2* Circular sample represented by unit 
masses and their center C o£ gravity. 

V — X cos «i, W — X sin at (4) 
The length of the resultant vector is 

R=(V^ + m)^/^ (5) 

Then the length r of the mean vector is 
simply 

r = -i-R (6) 

If the population center of gravity falls 
into the origin, then no single preferred direc- 
tion exists. This may be the case where all 
directions are equally likely, which is referred 
to as uniform distribution. Another possibility 
occurs with certain multi-modal distributions, 
for instance, when two opposite directiofts are 
equally probable. 

We should also consider the influence of 
chance fluctuations. A sample may have a 
non-zero mean vector even if the underlying 
theoretical distribution has its center of grav- 
ity in the origin. It is therefore desirable to 
test whether r differs from zero significantly. 
Such tests will be presented in the following 
five sections. 

Assume now that the length of the mean 
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vector is significantly different from zero and 
that the sample points are concentrated 
around a single preferred direction. Then we 
have good reason to introduce a mean direc- 
tion, which is defined to be the direction of 
the mean vector. We denote its angle by 
(Greek theta) and call it the mean angle 
(fig. 2) . This angle is calculated by solving 
the equations 


x — r cos ■&, y ==:r sin ^ 

(7) 

F == ii cos W — R sin ^ 

(8) 

RAYLEIGH TEST AND 


MODIFICATION 



The Rayleigh test is presented in reference 
1, p. 28. For reasons that are explained there, 
the Rayleigh test should be used only in the 
unimodal case. The null hypothesis states that 
the theoretical distribution is uniform. As a 
test statistic we may use the length r of the 
mean vector. When r exceeds a certain 
critical value, the null hypothesis is rejected. 

For the presentation of a table of critical 
values it is better to consider the test statistic 

z—nr^ (9) 

instead of r itself. A short table of critical 
values adapted from work by Greenwood and 
Durand was published in reference 1. Now 
it is possible to enlarge this table (table 1) 
considerably by adapting tables published in 
reference 2 and by using an unpublished table 
kindly submitted by W. T. Keeton, Cornell 
University. 

The Rayleigh test is most powerful if the 
alternative to the uniform distribution is a 
circular normal distribution (ref, 3) , 

Recently, a modification of the Ray- 
leigh test has proved to be most useful.^ 

^ I am indebted to Keeton for drawing my at- 
tention to the V test. 


It may occur that a particular direction is 
expected to be the preferred direction in ad- 
vance of the experiment. For instance, when 
pigeons are released at a test site, the home 
direction is known in advance. The null 
hypothesis that we are going to test is random- 
ne^, which means that the angles of the 
sample are independent observatioxis from a 
uniform circular distribution. For a t^t of 
the null hypothesis, it would be a loss of 
information when the knowledge of a pre- 
dicted direction were abandoned. Indeed, by 
using this direction we obtain a more pow- 
erful test. 

Let 00 be the angle of the predicted direc- 
tion, and let R, i? be defined as above. When 
making the predicted directions to a new 
X-axis, then the X-component of R is 

F' = R cos ^ (10) 

as shown in figure 3, By means of formula 
(8) we also get 

cos $o + W sin $0 (11) 

which is more practical for numerical cal- 
culations. 



FIGURE 3. The F' component of R as used In 
the V test. 
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Table 1,— Critical Values of the Test Statistic z of the Rayleigh Test 


n 

II 

O 

5 % 

2 . 5 % 

1 % 

0 . 1 % 

4 

2.36 

2.87 

3.28 

3.69 


5 

2.29 

2.84 

3.34 

3.86 

4.91 

6 

2.29 

2.86 

3.40 

4.08 

5.30 

7 

2.29 

2.88 

3.45 

4.16 

5.56 

8 

2.29 

2,90 

3.48 

4.20 

5.74 

9 

2.29 

2.91 

3.50 

4.25 

5.88 

10 

2.29 

2.92 

3.52 

4.29 

6.00 

11 

2.29 

2.93 

3.54 

4.32 

6.08 

12 

2.29 

2.93 

3.55 

4.34 

6.16 ■ 

13 

2.29 

2.94 

3.57 

4.36 

6.22 

14 

2.29 

2.94 

3.58 

4.38 

6.27 

15 

2.29 

2.95 

3.59 

4.40 

6.32 

16 

2 v 29 

2.95 

3.59 

4.41 

6.35 

17 

* 2.29 

2.95 


4.42 

6.39 

18 

' 2.29 

2.95 


4.43 

6.42 

19 

2.29 

2.96 

3.61 

4.44 

6.45 

20 

2.30 

2.96 


4.45 

6.47 

21 

2.30 

2.96 ‘ 

3.61 

4.46 

6.49 

22 

2.30 

2.96 

3.61 

4.46 

6.51 

23 

2.30 

2.96 


4.47 

6.53 

24 

2.30 

2.96 

3.63 

4.48 

6.54 

25 

2.30 

2.97 

3.63 

4.48 

6.56 

30 

2.30 

2.97 

3.64 

4.49 

6.62 

35 

' 2.30 

2.97 

3.64 

4.51 

6.66 

40 

2.30 

2.98 

3.65 

4.52 

6.69 

45 

2.30 

2.98 

3.65 

4.53 

6.72 

50 

2.30 

2.98 

3.66 

4.54 

6.74 

00 

2.31 


4.06 

4.61 

6.90 


In an experiment of homing, F' denotes 
the so-called homeward component (ref. 1, 
p. 19), if do is the angle of the homeward 
direction. 

The basic idea of the test is to consider 
the size of F'. If F' is small there is no evi- 
dence that the animals are oriented in the 
predicted direction with* angle do- If F' is 
relatively large, however, there must be some 
concentration of the directions around the 
predicted bearing. The larger the component 


F' is, the better chance there is of rejecting 
the null hypothesis of randomness. There- 
fore, we may choose F' as our test statistic. 
For this reason the test is called the F test 
(ref. 4). 

The V test leads to significance only if 
there is sufficient clustering around the pre-- 
dieted direction. In contrast the Rayleigh test 
is less powerful in this case but remains pow~ 
erful for clustering on any part of the circle, 

For preparing a numerical table of criti- 
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Table 2.-^Critical Values of the Test Statistic u of the Modified Rayleigh Test (V Test) 


n 

F - 0.10 

0.05 

0.01 

0.005 

0.001 

0.0001 

5 

M = 1.3051 

1.6524 

2.2505 

2.4459 

2.7938 

3,0825 

6 

1.3009 

1.6509 

2.2640 

2.4695 

2.8502 

3.2114 

7 

1.2980 

1 .6499 

2.2734 

2.4858 

2.8886 

3.2970 

8 

1.2958 

1.6492 

2.2803 

2.4978 

2.9164 

3,3578 

9 

1.2942 

1.6486 

2.2856 

2.5070 

2.9375 

3.4034 

10 

1.2929 

1.6482 

2.2899 

2.5143 

2.9540 

3.4387 

11 

1.2918 

1.6479 

2.2933 

2.5201 

2,9672 

3.4669 

12 

1.2909 

1.6476 

2.2961 

2.5250 

2.9782 

3.4899 

13 

1.2902 

1.6474 

2.2985 

2.5290 

2.9873 

3.5091 

14 

1.2895 

1 .6472 

2.3006 

2.5325 

2.9950 

3.5253 

15 

1.2890 

1.6470 

2.3023 

2.5355 

3.0017 

3.5392 

16 

1.2885 

1.6469 

2.3039 

2.5381 

3.0075 

3.5512 

17 

1.2881 

1.6467 

2.3052 

2.5404 

3.0126 

3.5617 

18 

1.2877 

1.6466 

2.3064 

2.5424 

3.0171 

3.5710 

19 

1.2874 

1.6465 

2.3075 

2.5442 

3 . 0211 . 

3.5792 

20 

1.2871 

1.6464 

2.3085 

2.5458 

3.0247 

3.5866 

21 

1.2868 

1.6464 

2.3093 

2.5473 

3.0279 

' 3.5932 ' 

22 

1.2866 

1.6463 

2.3101 

2.5486 

3.0308 

3.5992 

23 

1.2864 

1.6462 

2.3108 

2.5498 

3.0335 

3.6047 

24 ; 

1.2862 

1.6462 

2.3115 

2.5509 

3.0359 

3.6096 

25 1 

1.2860 

1.6461 

2.3121 

2.5519 

3.0382 

3,6142 

26 

1.2858 

1.6461 

2.3127 

2,5529 

3.0402 

3.6184 

27 

1.2856 

1.6460 

2.3132 

2.5538 

3.0421 

3.6223 

28 

1.2855 

1.6460 

2.3136 

2.5546 

3.0439 

3.6258 

29 

1.2853 

1.6459 

2.3141 

2.5553 

3.0455 

3,6292 

30 

1.2852 

1.6459 

2.3145 

2.5560 

3.0471 

3.6323 

40 

1.2843 

1.6456 

2.3175 

2.5610 

3.0580 

3,6545 

50 

1 .2837 

1 .6455 

2.3193 

2.5640 

3,0646 

3.6677 

60 

1.2834 

1.6454 

2.3205 

2.5660 

3.0689 

3.6764 

70 

1.2831 

1 .6453 

2.3213 

2.5674 

3.0720 

3.6826 

100 

1 .2826 

1.6452 

2.3228 

2.5699 

3.0775 

3.6936 

500 

1 .2818 

1.6449 

2.3256 

2.5747 

3.0877 

3.7140 

1000 

1 .2817 

1 .6449 

2.3260 

2.5752 

3.0890 

3.7165 


cal values, it is more practical to use the 
related test statistic 

u={2/n)^ F' (12) 

A chart of critical values published in refer- 
ence 4, p. 234, satisfies most practical 
purposes. The more accurate numerical table 
(table 2) was submitted by Keeton. 

For applications of the V test see refer- 
ences 5 and 6. The following examples were 
also suggested by Keeton: 


Example L Assume a random sample of 
directions as shown in figure 4. Despite an 
apparent clustering to the left, the Rayleigh 
test does not lead to significance (z =: 2.86, 
P > 0.5) . But when we learn that the home 
direction has azimuth — the pros- 

pect improves. Let us apply the V test. The 
azimuths oti^ listed in increasing order of 
magnitude, are 0°, 175®, 195®, 225®, 240®, 
240®, 260®, 295®, 330®, 340®, and 345®. 
With n = 11 we get, from formulas (4), (5) , 
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N 



FIGURE 4. Home direction is used as predicted 
direction in application o£ V test. 



FIGURE 5. Long experienced mean direction 
serves as predicted direction in application of V 
test. 

and (8), F= + 0.3514, = — 5.6026, 

/^ == 5.614, and t9^=274®. Then it follows 
from formulas (11) and (12) that F' = 
“(- 5.589 and u = 2.38. For n == 11, table 2 
reveals that P < 0.01. Hence the hypothesis 
of randomness can be rejected, whereas the 
Rayleigh test was too weak in our case. 

Example 2. At a test site near Castor 
Hill, N. Y., it was observed in a long series 
of releases that the mean bearing of homing 
pigeons always deviated by roughly the same 
angle from the homeward direction. Figure 5 


depicts the situation. The sample shown was 
collected under special experimental condi- 
tions that tend to increase the scatter of the 
bearings. The azimuths are 140®, 190®, 220®, 
230®, 255®, 270®, 300°, 330®, 330®, and 
350®. If we are interested in rejecting ran- 
domness, we may accept the long-experienced 
mean direction with azimuth Qq — 267.4® as 
a predicted direction. The test statistic of the 
V test is u — 2.21. Thus the V test leads to 
significance with a P value of less than 5 
percent, whereas the Rayleigh test would fail 
again since P turns out to be nearly 10 per- 
cent (z==2.45). 

TEST BY HODGES AND AJNE 

Hodges (ref. 7) proposed a bivariate sign 
test that later turned out to be quite useful 
in the analysis of directions. 

Assume that two quantities are measured 
on each individual or experimental unit. De- 
note the two measurements by x and y and 
assume that the measurements are repeated 
n times. Thus we get the pairs 

If the experiment is performed with the 
same individuals under different conditions, 
we get another sample consisting of n pairs, 
say 

{x{,yl) i ==l, 2, 

The null hypothesis states that the two 
samples were taken from the same bivariate 
distribution. To test this hypothesis we form 
the differences 

Xi — Xi, y{ — yiy i = 1, 2, • • • , n 

and plot these differences in a rectangular 
Xj );-diagram. Thus we get n vectors with 
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FIGURE 6. Hodges’ test for paired bivariate 
samples. 


components x{ ~ Xi and y{ — respective- 
ly (fig. 6) . If the null hypothesis is true, we 
would observe vectors pointing in all direc- 
tions. If, on the other hand, the null hypoth^ 
esis is wrong, the vectors should point to one 
sector of the x, y-diagram more frequently 
than to the rest of the plane. 

To get a suitable test statistic we may 
draw a straight line I through the origin and 
count how many of the vectors point to one 
or the other side of I, Now we rotate I 
around the origin until we minimize the 
number of vectors on one side. This minimum 
is the test statistic. We denote it by K. If K 
is small enough compared with the sample 
size n, the null hypothesis can be rejected. 
Critical values for K can be found in table 3, 
which is adapted from tables in references 3 
and 7. Hodges’ test may be interpreted as an 
extension of the sign test to bivariate data. 

Example 3. In figure 6 the sample size is 
n = 16, and the minimum number of vec- 
tors on one side of a straight line Hs = 2. 
The null hypothesis states that the measure- 


ments {xi, yi) and (x/, y/) belong to the 
same distribution. 

Table 3 yields F = 0.044. Hence we can 
reject the null hypothesis at a 5-percent level 
of significance. 

Example 4. In the paper in these pro- 
ceedings entitled “Satellite and Ground Radio 
Tracking of Elk,” F. C. Craighead et al. 
study the deviations between the actual loca- 
tion of an individual elk and the points 
tracked by a satellite. The errors of location 
are recorded as vectors in a horizontal plane 
in their figure 5. 

Do the vectors indicate any preferred di- 
rections? Applying Hodges’ test, we obtain 
n == 17 and K = 6, Table 3 yields a P-value 
greater than 0.500 or 50 percent. Hence we 
have no reason to assume a preferred direc- 
tion.2 

Notice, however, that east-west errors 
seem to be larger than north-south errors. To 
test significance of this kind we would need 
another test. Hodges’ test is not sensitive to 
changes of variance. 

Example 5. When tracking a single ani- 
mal, a question frequently arises: Is the 
movement of the animal oriented or is it 
random? Consider figure 7 where the track 
of an animal is recorded. At a glance it ap- 
pears that the animal is generally headed 
north and that it only occasionally deviates 
strongly from his course. 

In the light of probability theory, the 
track may be considered as a realization of a 
stochastic process. It is very likely that con- 
secutive sections of the track are dependent 
on each other. We want to show that the 
animal reorientates himself toward the north 
over and over again. 

Present-day statistics is not able to 


® I am indebted to C. E. Cote, Goddard Space 
Flight Center, Greenbelt, Md., for the permission 
to use his data. 
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Table Z — Critical Values of T est Statistics K in Relation to the Significance Level P and the 
Sample Size n 

n K =0 1 2 3 4 5 6 7 8 


8 P =^.063 .375 

9 .035 .246 

10 .020 .156 


11 

.011 

.097 

.376 







12 

.006 

.059 

.258 







13 

.003 

.035 

.171 

.489 



> .500 



14 

.002 

.021 

.111 

.355 






15 

.001 

.012 

.070 

.250 






16 


.007 

.044 

.171 

.444 





17 


.004 

.027 

.114 

.327 





18 


.002 

.016 

.075 

.233 





19 


.001 

.010 

.048 

.163 

.399 




20 



.006 

.030 

.111 

.296 




21 



.003 

.019 

.074 

.213 

.466 



22 



.002 

.012 

.049 

.151 

.356 



23 



.001 

.007 

.032 

.104 

.265 



24 




.004 

.020 

.071 

.193 

.413 


25 

A 

*o 

o 



.003 

.013 

.048 

.137 

.315 


26 




.002 

.008 

.031 

.096 

.235 

.466 

27 





.005 

.020 

.066 

.172 

.364 

28 





.003 

.013 

.045 

.124 

.278 

29 





.002 

.008 

.030 

.087 

.208 

30 





.001 

.005 

.020 

.061 

.153 

35 







.002 

.008 

.026 

40 









.003 


prove that such a track is oriented. Nor is our example^ we assume that a new decision is 

statistics able to deal with random walks, made at intervals smalto than an l^r.) Now 

However, the following modest approach is we form the vectors PoPi, . 

possible: Let us play the role of a devil’s Under the null hypothesis these vectors have 

advocate and assume that from time to time random directions. We plot these vectors 

the animal chooses new directions ‘^at ran- with a common base 0 and draw a line t 

dom.” The directions would be taken from a that minimizes the number of vectors on one 

uniform circular distribution. We call this side of Z. Applying Hodges’ test to the vectors 

assumption the null hypothesis. We consider of our example, we get n = 13, K — I, 

the points Pi, Ps, • " * reached by the animal Table 3 yields P ==s .035. Hence we are able to 

at time 1 hour, 2 hours, • * • j respectively, reject the null hypothesis of randomness at a 

after release at Pq. (The time interval de- level of P = 5 percent, 

pends on the frequency of new decisions. In Ajne (ref. 3) proposes a circular one- 
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FIGURE 7. Track of an animal and a test for 
randomness. 

sample test. The problem is the same as for 
the Rayleigh test discussed previously. Con- 
sider a sample from a unimodal circular dis- 
tribution. Can a preferred direction be 
established statistically? The null hypothesis 
is a uniform distribution. 

The test procedure is very simple once a 
circular plot is provided. We draw a straight 
line I through the center of the circle and 
count the number of sample points on each 
side of L Then we rotate I and minimize the 
number of sample points on one side of Z. 
This minimum number is denoted by K and 
used as a test statistic. If K is small relative 
to the sample size n, the null hypothesis of a 
uniform distribution is rejected. 

We slightly modified the original formu- 
lation of Hodges’ and Ajne’s tests. Thus it 
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becomes obvious that Ajne’s test is a special 
case of Hodges’ test. This connection was dis- 
covered by Bhattacharyya and Johnson (ref. 
8) . Ajne’s test is powerful if the alternative 
to the uniform distribution is a rather narrow 
unimodal distribution (high degree of con- 
centration around the mean direction) . 

Example 6. Homing pigeons were released 
singly. They disappeared at directions meas- 
ured by the following azimuths (arranged in 
increasing order) : 115°, 120®, 120®, 130®, 
135®, 140®, 150®, 150®, 150®, 165®, 185®, 
210®, 235®, 270®, and 345®. Sample size is 
n = 15. The line I can be drawn in such a 
way that the minimum number of sample 
points on one side of Z is = 1 (fig. 8) . 
From table 3 we get P = .012. Hence the 
null hypothesis of a uniform distribution can 
be rejected at a 2-percent level of significance. 

AJNWS SECOND TEST 

Ajne (ref. 3) proposes a further method 
for testing the null hypothesis of a uniform 
distribution. Consider again the straight line 
Z in figure 8. If the theoretical distribution 
were uniform, we would expect an equal 



FIGURE 8. Ajne’s test applied to an orientation 
problem. 
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Table 4:,— Critical Values for Ajne^s Test Statistic A in Degrees 


n 

P=10% 

5% 

2.5% 

1% 

0.5% 

5 

A = 185 

227 

262 

301 . 

324 

6 

184 . 

227 

268 

314 

343 

7 

184 

228 

269 

318 

354 

8 

184 

229 

271 

322 

359 

9 

184 

230 

272 

326 

364 

10 

185 

231 

274 

329 

369 

11 

185 

231 

275 

332 

372 

12 

185 

232 

277 

334 

375 

13 

185 

232 

277 

336 

377 

14 

185 

232 

278 

337 

379 

15 

185 

232 

278 

338 

381 

16 

185 

233 j 

279 

339 

383 

17 

185 

233 ^ 

279 

340 

385 

18 S 

185 

233 

280 

341 

386 

19 

185 

233 

280 

341 

386 

20 

185 

233 

281 

342 

387 

30 

186 

234 

283 

346 

393 

40 

186 

235 

284 

347 

395 

50 

186 

235 

284 

349 

397 

100 

186 

235 

285 

351 

400 

200 

186 

236 

286 

352 

402 

00 

00 

236 

287 

354 

404 


number of sample points on each side of h 
that is, n/2 sample points. A small deviation 
from n/2 may be due to chance fluctuation. 
A large deviation, however, indicates that 
the uniform distribution is not the proper 
hypothesis. When we rotate I around the 
center, the deviation from n/2 varies. Ajne 
chooses a suitably defined average of this 
deviation as a test statistic. 

A rather lengthy mathematical analysis 
leads finally to the following procedure (ref. 
9). Let «!, • • • 3 «n be the observed 

angles in degrees and assume that they are 
arranged in ascending order. Thus 


Calculate the differences 

mi2 = a2 — «i, mi3 = — ai, 

m23 = «3 — • ‘ • 

If one of them exceeds 18G°, take the comple- 
mentary value by subtracting it from 360®. 
For instance, if aio — 205®, take 

mi,io = 360® — 205® = 155®. The differ- 
ences may be arranged as follows: 

mi2 mi3 mi4 mis ’ ’ ' rn^^n 

W23 ?W24 m25 • • • rri2,n 
ms5 : 

W245 : 


ai 02 «3 ^ ( 13) 


(14) 
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Then form the sum of all these n(n — 1) /2 
differences, that is, 

Z= i mij (15) 

t— 1 

With this sum the test statistic can be writ- 
ten in the form 

A = n-90°- —Z. (16) 

n 

The null hypothesis of a uniform distribu- 
tion is rejected if A exceeds a certain critical 
value. A table of critical values was published 
by Stephens (ref. 9) ; our table 4 is an adapta- 
tion of it. 

Ajne’s second test is especially powerful 
if the alternative to the uniform distribution 
is a unimodal distribution with a large angular 
deviation (low degree of concentration) • 

Example 7. We apply Ajne’s test to the 
same data as in example 6. Here a\ = 1 15 ® , 
0(2 = 1 20 ° , * • * , ais = 345 ® . The differences 
in degrees arranged in the scheme (equation 
14) are 


5 

15 

20 

25 

35 

35 

35 

0 

10 

15 

20 

30 

30 

30 


10 

15 

20 

30 

30 

30 



5 

10 

20 

20 

20 




5 

15 

15 

15 





10 

10 

10 


0 0 • • • 

0 • * • 

Notice thata 55 — Ck:i = 345^ — 115^ = 230^^, 
but mi,i 5 = 360® — 230® = 130®. Formula 
(15) yields Z = 6660®. Hence, we get from 
formula (16) 

4 = 15 (90 ° ) — 6660° = 462° 


This value is larger than any critical value in 
table 4. Therefore, we can reject the null 
hypothesis of a uniform distribution at a 0.5- 
percent level of significance. 

TEST BY LAUBSCHER AND RUDOLPH 

Sometimes tests are required that are not 
necessarily powerful but are quick to apply. 
When testing uniformity versus a unimodal 
alternative on the circle, such a test could 
be quite useful. This was proposed by Laub- 
scher and Rudolph (ref. 10) and also in- 
vestigated by Rao (ref. 11). 

Let i? be the length of the smallest arc on 
the circle that contains all sample points. R 
may be called the range of the sample. A 
sufficiently small value of R indicates that the 
sample was not taken from a uniform dis- 
tribution. Therefore, R can be chosen as a 
test statistic. If R is below a certain critical 
value, the null hypothesis of a uniform dis- 
tribution can be rejected. A table of critical 
values follows (table 5) which is part of a 
table published in Laubseher and Rudolph 
(ref. 10). 

Example 8. Six different migrating birds 
kept in cages moved toward average direc- 
tions with the following azimuths : 122®, 93®, 
158®, 67®, 85®, and 145®. Is the concentra- 
tion of directions significant? 

All six angles are between 67® and 158®. 
Hence the range is 12 = 158® — 67® = 91®. 
The critical value of R for a significance 
level of 1 percent is 100.2® (table 5). The 
sample range is below the critical value. 
Hence the null hypothesis of a uniform dis- 
tribution can be rejected. The concentration 
of directions is indeed significant. 

RAO’S TEST 

The previous sections have dealt with the 
same topic— testing whether a circular sam- 
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Table 5,— Critical Values n Degrees of the Range K of a Circular Sample 


n 


0.01 

0.025 

0.05 

0.10 

4 

R-38.8 

48.9 

66,3 

83.6 

105.3 

5 


76.1 

95.7 

113.8 

135.4 

6 

87.2 

100.2 

120.3 

138.2 

158,7 

7 

107.6 

120.8 

140.8 

158.0 

177.3 

8 

125.5 

138.5 

157.9 

174.4 

192.5 

9 

141.1 

153.8 

172.5 

188.1 

205.1 

10 

154.7 

167.1 

185.0 

199.8 

215.8 

11 

166.8 

178.7 

195.9 

209.9 

225.0 

12 

177.4 

189.0 

205.4 

218.7 

233.0 

13 

187.0 

198.1 

213.8 

226.5 

240.0 

14 

195.5 

206.2 

221.3 

233.4 

246.2 

15 

203.2 

213.5 

228.0 

239.5 

251.7 

16 

210.2 

220.1 

234.0 

245.1 

256.7 

17 

216.6 

226.2 

239.5 

250.1 

261.2 

18 

222.4 

231.6 

244.5 

254.7 

265,3 

19 

227.7 

236.7 

249.0 

258.8 

269.1 

20 

232.7 

241.3 

253.2 

262.7 

272.5 

21 

237.2 

245.6 

257.1 

266.2 

275.6 

22 

241.4 

249.5 

260.7 

269.4 

278.6 

23 

245.4 

253.2 

264.0 

272.5 

281.3 

24 

249.0 

256.7 

267.1 

275.3 

283.8 

25 

252.5 

259.9 

270.0 

277.9 

286.1 

26 

255.7 

262.9 

272.7 

280.4 

288.3 

27 

258.7 

265.7 

275.2 

282.7 

290.4 

28 

261.5 

268.3 

277.6 

284.8 

292.3 

29 

264.2 

270.8 

279.8 

286.9 

294.1 

30 

266.7 

273.2 

281.9 

288.8 

295.9 


pie is taken from a uniform distribution. All 
of these tests are restricted to the situation 
where the alternative is a unimodal distri- 
bution. 

However, multimodal distributions occur 
frequently in biology.^ They have two or 
more preferred directions. When we apply 
one of the aforementioned tests to such data, 
the test loses most of its power; that is, it 


*We will omit here the special treatment of 
multimodal distributions with strict symmetry. 


frequently fails to reject the null hypothesis 
of a^ uniform distribution even though the 
existence of modes is strongly suggested. 

A test that is powerful in unimodal as 
well as multimodal situations was recently 
proposed by Rao (ref. 11). Its application is 
easy. Assume that the angles of a random 
sample are arranged in increasing order: 

ai ^ a2 — «3 — • • • an 
Then we calculate the length of all n arcs 
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Table 6,-^Critical Values for the Test 
Statistic U in Degrees 


n 

P=0.01 

0.05 

0.10 

4 

U = 221.0 

186.5 

171.7 

5 

212.0 

183.6 

168.8 

6 

206.0 

180.7 

166.3 

7 

202.7 

177.8 

164.9 

8 

198.4 

175.7 

163.4 

9 

195.1 

173.5 

162.4 

10 

192.2 

172.1 

161.3 

11 

189.7 

170.3 

160.2 

12 

187.6 

169.2 

159.2 

13 

185.8 

167 .8 

158.4 

14 

184.0 

166.7 

157.7 

15 

182.2 

165.6 

157.0 

16 

180,7 

164.9 

156.6 

17 

179.6 

164.2 1 

155.9 

18 

178.2 

163.1 1 

155.2 

19 

177.1 

162.4 

154.8 

20 

176.0 

161.6 

154.4 

25 

171.9 

158.9 

152.7 

30 

168.8 

156.7 

151.4 

35 

166.4 

155.0 

150.3 

40 

164.4 

153.6 

149.5 

45 

162.7 

152.4 

148.7 

50 

161.2 

151.4 

148.1 

100 

152.8 

146.8 

143.7 

200 

146.8 

142.6 

140.4 


between consecutive sample points on the 
circle. We denote these arc lengths by Tt. 


— ~ au 

T2 = 0:3 — CK23 


Tn— 1— Otn 

r„ = 360° + ai — (17) 

Notice that we always have 

2 ^,^360^ (18) 



FIGURE 9* Toting a bimodal sample using 
Rao^s test. 


we expect that the differ only slightly 
from each other. They fluctuate around their 
mean value 360® /n. On the other hand> if 
the Ti differ sufficiently from 360® /n^ this 
indicates that the theoretical distribution is 
not uniform but rather unimodal or multi- 
modal. Rao introduces the sum of deviations 
as test statistic ; more specifically 



n 

2 

=: 1 



(19) 


If U exceeds a certain critical value, the 
null hypothesis of a uniform distribution is 
rejected. The following (table 6) is a table 
of critical values for the U statistic. It is 
partly based on a table in reference 1 1 . 

Example 9. Homing pigeons were re- 
leased singly in the Toggenburg Valley under 
subalpine conditions.^ The birds did not ad- 
just quickly to the homing direction but pre- 
ferred to fly in the axjs of the valley (fig. 9) . 
The vanishing points are given by the angles 
arranged in ascending order: 20®, 135®, 
145®, 165®, 170®, 200®, 300®, 325®, 335®, 
350®, 350°, 350®, and 355®. Sample size is 
n=13. 


, , * I am indebted to G. Wagner, Bern, Switzerland, 

If the theoretical distribution is uniform, for the permission to use his data. 
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Table 7.- — Critical Values of the Test Statistic K = y/n (i?+ + D~) 


n 

P =10% 

5% 

1% 

0.1% 

4 

K = 1.428 

1 .536 

1.728 

1.874 

5 

1.458 

1.565 

1.764 

1.970 

6 

1.473 

1.584 

1,793 

2.018 

7 

1.484 

1.598 

1.815 

2.051 

8 

1.493 

1 .609 

1.830 

2.076 

9 

1.501 

1.618 

1.842 

2.096 

10 

1.507 

1 .625 

1.853 

2.112 

11 

1.512 

1.632 

1.862 

2.126 

12 

1.518 

1.637 

1 .869 

2.140 

14 

1.527 

1.647 

1.882 

2.155 

16 

1.534 

1.656 

1.892 

2.168 

18 

1.540 

1.662 

1.901 


20 

1.546 

1.665 

1.908 


30 

1.562 

1 .684 

1.930 


40 1 

1 - 571 

1.695 

1 .941 


50 

1.577 

1.701 

1-949 


60 

1.582 

1.705 

1.955 


70 

1.585 

1 .708 

1.959 


80 

1.588 j 

1.711 

1.962 


100 

1.590 

1.716 

1 , 967 


oo 

1.620 

1 .747 


2. 303 


Does a uniform distribution fit these 
data? To test this possibility we apply Rao’s 
U statistic. From formulas (17) Xve obtain 

Ti=115° T5= 30° r9==15° 

T 2 = 10 ° r6 = ioo° rio= o° 

T3= 20° r7== 25° rii= 0° 

T4= 5° T8= 10° 7x2= 5° 

7x3 = 25° 

As a check of computation we get X7t 
360°. 

Since 360°/n = 27.7°, we obtain from 
formula (19) 

C7==-|-(87.3 4- 17.7 +7.7 + 22.7 + 2.3 
+ 72.3 + 2.7 + 17.7 + 12.7 + 27.7 
+ 27.7 + 22.7 + 2.7) = 162° 


For P == 0.05 we get from table 6 the 
critical value 167.8^; our U value from the 
sample is somewhat below this. Thus we 
were not able to establish significance at a 
5-percent level. However, a slightly higher 
sample size^^ay n = 15, would probably have 
yielded the desired result. 

KUIPER’S GOODNESS^OF-^FIT TEST 

In reference 1, Kuiper’s test was described 
in detail. The test could serve the same pur- 
pose as the previous tests, that is, to decide 
if a circular distribution is uniform. 

Stephens (ref. 12, table 1) slightly cor- 
rected the table of critical values for the test 
statistic and enlarged it considerably. Table 
7 is based on Stephens’ findings. 
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Kuiper’s test can be used in the unimodal 
and the multimodal case. A test related to 
Kuiper’s test is Watson’s W test. It is de- 
scribed in reference 1 (p. 27) . In samples 
occurring in biology, Watson’s test is at least 
as powerful as Kuiper’s test. 

Example 10, We apply Kuiper’s test to 
the data of example 9. Following the test 
procedure as explained in reference 1, p, 26, 
we obtain figure 10. From this figure we read 

D+ == .023 and D” = .369. Hence 
K = + D- ) = *023 + .369) 

= 1.414 

From table 7 we obtain the critical value 
K = 1.642 at a 5-percent level of significance. 

Our K value from the sample is smaller. 
Therefore, we cannot claim a significant 
deviation from the uniform distribution. The 
result is consistent with our findings based oh 
Rao’s test. 


CONFIDENCE INTERVALS FOR 
THE MEAN ANGLE 

In the case of a unimodal distribution, it 
is convenient to work under the hypothesis 
of circular normal distribution introduced by 
von Mises in 1918. The probability density 
function is 

f{a) = (20) 

where 6 denotes the mean angle, K the pa- 
rameter of concentration, and C a numerical 
constant depending on k. More details, charts, 
and tables are given in reference 1. 

If we are given a sample of n independent 
angles 0 Jx,a 2 , " * * it is often demanded 
to estimate the parameters 6 and k. As a 
point estimate of 0, we take the mean angle 
calculated from the sample (formula 7) . 

In addition, the research worker wants to 
know the reliability of this estimate. The 



JFIGURE 10. Application of Kuiper’s test to a circular sample. 
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n S= 200 100 50 40 30 20 12 10 8 



FIGURE 11. Chart for determining a confidence interval of the mean angle with a 95 percent 
confidence co^dent. 


question is then to find an interval around ^ 
in which the unknown parameter 6 falls with 
a preassigned probability called the con- 
fidence coefficient. 

In reference 1 a confidence interval of 
the form 

^±8 (21) 

was presented. However, the method to de- 
termine the deviation 8 is somewhat sophis- 


ticated. It was therefore a good idea (ref. 
13) to search for a direct graphical approach. 
Brown and Mewaldt present a chart from 
which the deviaton 8 can be read as a func- 
tion of the sample size n and the length r of 
the mean vector. Their chart was prepared 
for only Q = 99 percent and is fairly sketchy. 
It seemed worthwhile to follow their idea 
and to draw accurate charts. Figures 1 1 and 
12 yield 8 for confidence coefficients Q = 95 
percent and Q = 99 percent, respectively. 
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FIGURE 12. Chart for determining a confidence interval of the mean angle with a 99 percent 
confidence coefficient. 


Example IL In example 6 the null hy- 
pothesis of a uniform distribution was re- 
jected. It is reasonable to assume that the 
sample was taken from a circular normal 
distribution with unknown mean angle d 
and unknown parameter k of concentration. 

To estimate B we first determine the cen- 
ter of gravity following the instructions given 
in the second section of this paper. We get 

2 cos = — 8.5723, 


8.5723) = — .5715, 

% sin =: -f- 3.8478, 

y = ^(3.8478) = + ,2565, 

r = (-f2 + y2) = .6264, 

co$d = — .9123, 

r 

sin # ^ = 4- .4095, 

155.8° 

Now we choose Q = percent as con- 
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0 j .2 .3 .4 .5 .6 .7 .8 .9 1.0 


length of mean ve c t or r 

flGURE 13. Chart for determiniiig a confidence interval of the parameter of concentration with 
a 90 percent confidence coefficient. 

fidence coeflBicient. In figure 1 1 we find lines 
for sample sizes n == 14 and 16. Our sample 
size is n == 15 so that we have to interpolate. 

At r = .6264 we find 8 — 30.6'^ for n = 16 
and 8 = 33.4° for n = 14. Hence 8 = 32.0° 
fits our purpose. Using formula (21) our con- 
fidence interval for 6 turns out to be 

155.8° ± 32.0° 

CONFIDENCE INTERVALS FOR 
PARAMETER OF CONCENTRATION 

As in the previous section, we assume 


that the theoretical distribution is a circular 
normal distribution with probability density 
function (20). Our purpose is to estimate 
the parameter /c. 

Values of k range from 0 to oo . The par- 
ticular value K = 0 indicates that the dis- 
tribution is uniform. The higher the value 
of Ky the stronger is the concentration about 
the mean direction. Usually k has to be 
estimated when the mean angle $ is unknown. 
The estimation is based on the length r of 
the mean vector, which was defined in 
formula (3). For a point estimate of k, we 
use table B in the appendix of reference 1. 
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FIGURE 14. Chart for determining a confidence interval of the parameter of concentration with 
a 98 percent confidence coefficient. 


There r ranges from 0 to 1 in steps of 0.01. 
The estimate of k can he read in the third 
column. 

Stephens (ref. 14) studied the theoretical 
distribution of r for a variety of circular 
normal distributions where k ranges from 0 
to 5. Based on his table 2^ the following 
charts were drawn: Figures 13 and 14 give 
upper and lower confidence limits for k for 
various sample sizes n and for the confidence 
coefficients Q = 90 and 98 percent, respec- 
tively. The use of the charts is explained in 
example 12. 


We denote the lower limit by k% and the 
upper limit by K«. Thus the confidence in- 
terval is 

Ki<K< Ku (22) 

Notice that sometimes the lower confi- 
dence limit is zero so that there is no uncer- 
tainty on the lower side of the interval. In 
this case the confidence coefficient has to be 
raised approximately from 90 to 95 percent 
and from 98 to 99 percent, respectively. 

Example 12. We return to example 11. 
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There n = 15 and r == .6264. To get a point 
estimate of k, we use table B in the appendix 
of reference 1. For r=:.62 we find 
K — L60044. For r = ,6264 we get by linear 
interpolation k = 1 .629. 

For confidence limits of k, we choose 
(2 = 90 percent as our confidence coefficient. 
In figure 13 we find curves for n == 10 and 
n = 20. For our own sample with n =15 
we will have to interpolate. For the lower 
limit we read from the lower curves at 
r .6264 

71 = 10, Ki = 0,44, 

n = 20j Ki = 0.88, 

By linear interpolation for n==15 we 
obtain 

Ki = 0.66 

Similarly for the upper limit we read from 
the upper curves at r = .6264 

n — 10, Ku — 2.62, 

n = 20, Ku — 2.40. 


Linear interpolation for n = 15 leads to 



FIGURE 15, Testing two bivariate samples. One 
sample is represented by open circles, the other 
sample by filled circles. (Mardia’s test). 


Ku == 2.51 

Hence with probability 90 percent we con- 
clude that 

0.66 < K < 2.51 

Notice that the point estimate k= 1.629 
does not fall exactly into the center of this 
interval. 

TWO^S AMPLE TEST BY MARDI A, 
WATSON, AND WHEELER 

Wheeler and Watson (ref. 15) proposed 
a test procedure for comparing two circular 
samples. This test is nonparametric and pow- 
erful; for both reasons the test should be 
favored by research workers. 

Mardia (ref* 16) has shown that Wat- 
son’s and Wheeler’s test is a special case of 
a bivariate test proposed by Mardia (ref. 
17). Mardia also provided tables of critical 
values for the test statistic. Therefore, it 
seems to be appropriate to attribute the cir- 
cular test to all three authors and to name it 
the Mardia-Watson- Wheeler test. In this sec- 
tion we will first explain Mardia’s bivariate 
test and later consider the special case for 
circular samples* 

We consider a sample of bivariate obser- 
vations 

(xi,yi), i — 1,- 
and a second sample 

(•^Ay/)> i = 

We want to know whether the two samples 
belong to the same bivariate population. To 
explain the test procedure we plot the sample 
points (fig. 15). The sample points (xi, yi) 
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FIGURE 16. (a) Reduction of bivariate samples to circular samples, (b) Generation of equidis- 
tant sample points (Mardia-tVatson-Wheeler test). 

Now we rank the m n directions with 
numbers 1 , 2 , 3, • • V, m + n by starting at 
any reference direction and by rotating coun- 
terclockwise. Let 

be the ranks of the first sample^ and 
ri', • • • , 

the ranks of the second sample. In figure 16a 
we see that ri = 1 , r 2 = 2 , = r 4 = 4, 

Ts = 83 re = 9 and == 5, t 2 == 6 ^ r/ = 7y 
r/ = 10. We have tacitly assumed that there 
are no intersample ties. 

As a next step we introduce n -f- 
equally spaced dots on the circumference of 
the unit circle and mark them in the same 
order as the dots of the original samples {fig. 
16b) . If the dots of one sample are sufficiently 
separated from the dots of the other sample, 


of the first sample are plotted by filled circles 
and the sample points (x/, y/) of the second 
sample by open circles. The second sample is 
located slightly to the left and below the first 
sample. Is this shift due to chance only or 
were the samples drawn from different popu- 
lations? 

To find a statistical answer to this ques- 
tion we pool the two samples and calculate 
the coordinates of the common center C of 
gravity using the formulas 


%Xj + 

m n ’ 


( 23 ) 

m n 


The center C of gravity is marked in figure 
15. 

From C we draw vectors to all m n 
sample points and consider their directions. 
Mardia’s test is based on these directions 
only ; thus the test reduces the bivariate case 
to the circular case (fig. 16a). 
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we conclude a significant difference between 
the samples. To define a test statistic we 
concentrate on one of the two samples. The 
choice is irrelevant. Then we calculate the 
length R of the resultant vector using for- 
mulas (4) and (5). When choosing the first 
sample, the angles are 

r=l,2,---,m(24) 

n -f- m 

Hence 

m m 

V == % cos ^i,W = S sinjSi, 

(25) 

In figure 16b the angles )8i are plotted 
with the horizontal direction as zero line : 

^1 = 36°, iS2 = 72^ j 83 = lG8°, 

^4=144°, )85 = 288^ ^6 = 324° 

Hence 1.118, 1.539, i? — 1.902. 

This length R can be used as a test sta- 
tistic. If R is sufficiently large, the dots of the 
first sample are concentrated around a pre- 
ferred direction and thus more or less sepa- 
rated from the dots of the second sample. 
This argument is very similar to the pro- 
cedure of the Rayleigh test {see the third 
section) . Thus, in a certain sense, the two- 
sample test is reduced to a one-sample test, 
an idea which was successfully applied to 
other tests (ref. 18) . 

Instead of R itself, Mardia uses the fol- 
lowing function of R as a test statistic 

B = R^ (26) 

If B exceeds a certain critical value, we re- 
ject the null hypothesis that the samples were 
taken from the same population. 

The Mardia- Watson- Wheeler test can be 


safely used only if there are no intersample 
ties, that is, if no angle of the first sample 
coincides with an angle of the second sam- 
ple. When the test leads to significance, this 
does not necessarily imply that the two sam- 
ples differ in location. It could well occur 
that the main reason for significance is a dif- 
ference in dispersion. This property is typical 
for most nonparametric tests. For a discussion 
see reference 1 (p. 34, especially figure 23.1). 

Table 8 of B values is based on table 1 of 
reference 17. 

For N > 17, Mardia (ref. 17) gives the fol- 
lowing approximate distribution : 

The quantity 

mn ^ ' 

is approximately distributed as ^ with two 
degrees of freedom. Using a table of critical 

values, it is easy to find critical values for 
U. 

In our illustrative example of figures 15 
and 16 we got R= 1.902. Hence, B = R^ 
= 3.62. For n m =: N = 10 and n = 4, 
the critical value is 9.47 at a 5 percent level 
of significance. Our B value is smaller. Thus 
we cannot claim significance. 

For a biological example see Mardia (ref. 
16, p. 189). 

WATSON’S W TEST 

A nonparametric two-sample test attrib- 
uted to Cramer, von Mises, and Smirnov was 
adapted for the circular case by G. S. Wat- 
son. The test was explained and exemplified 
in reference 1 (pp. 35 and 36). The test sta- 
tistic is denoted by m where n and m are 
the two sample sizes. 

The lack of critical values of m has 
prohibited the use of this test in the begin- 
ning. Now tables of critical values are avail- 
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Table 8. — Critical Values of the Test Statistic B. (The Two Samples Are of Size m and n, 
m ^ N = m + n.) 

n P =0 . 001 0 .01 0 . 05 0 . 10 
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Table 9.— Critical Values for Watson^s U*n,m- (The Two Samples Are of Sizes n and m^n^m) 


n^m 

i »=10% 

5% 

1% 

0.1% 

n,m 

P =10% 

5% 

1% 

0.5% 

10,6 

.153 

. 182 

.243 

.266 

20, 4 

,151 

.178 

.231 

.251 

10, 7 

.153 

.183 

.246 

.270 

20,6 

.152 

.182 

.245 

.270 

10, 8 

.153 

.183 

.248 

.273 

20, 8 

.152 

.183 

.251 

,278 

10, 9 

, 153 

.184 

.249 

,275 

20, 12 

.153 

.185 

.256 

.285 

10, 10 

.153 

.184 

.250 

.277 

20, 16 

.153 

.185 

.258 

.289 

12, 6 

.153 

.182 

.244 

.268 

20, 20 

.153 

.185 

.259 

.290 

12,7 

.153 

.183 

.247 

.272 

25, 5 

.151 

.180 

.240 

.263 

12, 8 

.153 

.184 

.249 

.275 

25, 10 

.152 

.184 

.254 

.283 

12, 9 

.153 

.184 

.251 

.278 

25, 15 

.153 

.185 

.258 

.289 

12, 10 

.153 

.184 

.252 

,279 

25, 20 

.153 

,186 

.260 

,291 

12, 11 

.153 

.184 

.253 

.280 

25, 25 

.153 

,186 

.261 

.293 

12, 12 

.153 

.184 

.253 

.282 

30, 5 

,150 

.179 

.240 

.263 

14, 6 

.153 

.182 

.244 

.269 

30, 10 

.152 

.184 

.254 

.284 

14, 8 

.153 ' 

.184 

.250 

.276 

30, 15 ! 

.152 

.185 

.259 

,290 

14, 10 

.153 

.184 i 

.253 

.281 

30, 30 

,153 

.186 

.262 

.294 

14, 12 

.153 

.185 

.254 

.283 

40, 5 

.150 

,179 

.240 

.263 

14,14 

.153 

.185 

.255 

.285 

40, 10 

.152 

.184 

.254 

.284 

16,4 

.151 

.178 

.231 1 

.251 

40, 15 

, 152 

.185 

.259 

.290 

16,6 

.152 

.182 

.244 

.269 i 

40, 20 

.152 

.186 

,261 

.293 

16,8 

.153 

.183 

.250 

.277 

40, 40 

.152 

.186 

.263 

.296 

16, 10 

.153 

.184 

.253 

.282 

50, 5 

.150 

.178 

.239 

,263 

16, 12 

.153 

.185 

.255 

.284 

50, 10 

.151 

.183 

.254 

.284 

16,14 

. 153 

.185 

.256 

.286 

50, 15 

.152 

.185 

.259 

.290 

16, 16 

.153 

.185 

.257 

.287 

50, 20 

.152 

.185 

.261 

.293 






50, 25 

.152 

.186 

.262 

.295 






50, 50 

.152 

. 186 

.264 

.298 






00 , 00 

,1517 

.1869 : 

.2684 

.3035 


able. Burr (ref. 19, p. 1094) published a 
table for all pairs m, n ^ 4, m + n ^ 17, 
and P < 0.01. Later Stephens (ref. 20) , ex- 
tended the table for large values of n and 
m up to m — n = 50. 

Watson’s U\ m test is sensitive for all 
kinds of deviations between two populations. 
Table 9 is adapted from reference 20. 

Vn^mTEST 

Another nonparametric two-sample test, 
Kuiper’s modification of a test by Smirnov, is 
now known as the test. It is described in 
detail in reference 1 (p. 35) , Maag and 
Stephens (ref. 29) have provided tables for 


the critical values of the test statistic. They 
range from n = 3 to n = 10 and from 
m = 3 partly up to m = 20. 

MULTISAMPLE TESTS 

In orientation problems, sometimes more 
than two samples have to be compared with 
each other. Tests designed to discover differ- 
ences among several samples are often called 
tests of homogeneity, 

A test for comparing the mean directions 
of several circular samples was proposed by 
Watson and Williams (ref. 21) . It is assumed 
that q samples of sizes tzi, n 2 j . - . , Uq, re- 
spectively, are taken from q circular normal 
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. (or von Mises’) distributions with the same 
unknown parameter k of concentration. 

The null hypothesis states that the mean 
directions of the q normal distributions are 
the same. To test the null hypothesis we cal- 
culate jPor each sample the length Ri of the 
resultant vector using formulas (4) and (5). 
We also pool the q samples and calculate 
for all 

JV == n X *4“ ^2 "4" * ’ * *4" (28) 
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concentration. A nonparametric r^ultisample 
test was proposed by Maag (ref. 24). It is an 
ejctension of Watson’s IP test (see the previ- 
ous section). Unfortunately the test cannot 
yet be appjlied since a table of critical values 
for the test statistic is not available. 

A generalization of the Mardia-Watson- 
Wheeler test to the multisaniple case will 
soon be published. For this test critical values 
of the test statistic will be available (see refs. 
25 and 26). 


observations the length R of the resultant 
vector. Then the test statistic is 


lN~q)gR,-R) 


(29) 


If F exceeds a certain critical value, the null 
hypothesis of equal mean directions is re- 
jected. The test statistic F is approximately 
distributed as Fisher’s Fq^i, N—q with q — -I 
and N^q degrees of freedom. Critical 
values pf F can be found in any table of the 
F distribution, e.g. in the well-known table 
by Fisher and Yates (ref. 22) on pages 47, 
49, 51, 53, 55. There F is denoted by 

For the rationale of the test and for a 
biological example see reference 1 (section 
22 ). 

A modification of Watson’s and Williams’ 
test was recently proposed by Stephens (ref, 
23) . For the two-sample case, Stephens’ test 
is exact. In the multisample case, the accu- 
racy is improved by the introduction of a fac- 
tor G that depends on R/N, 

These tests are parametric; that is, they 
are based on strict assumptions on the under- 
lying distributions. They axe certainly power- 
ful, but the power may be outweighed by 
doubts about the basic assumptions. It is 
therefore desirable to apply nonparametric 
tests whenever we have no evidence of uni- 
modal distributions with the same degree of 


COSINOR METHOD : 

In the statistical analysis of biologic^^^ 
rhythms a method has been developed that is 
strongly connected to circular distributions. 
If the period T of a cyclic phenpinenon is 
known, it is natural to associate the quantity 
measured with a circle. As thp time increases 
by one period, we rotate once around t})k 
circle. As a rule, there will be a time instant 
when the quantity reaches a maximum. TKis 
time instance corresponds to the mean difec- 
tion of a unimodal circular distribution. . 

To be more specific, we consider the data 
in figure 17.® The interperitoneal tempera- 
ture of a rat is measured at various times 
during the 24 hours of a day, hot necessarily 
at regular time intervals. The data suggest 
the following model : 

Zi Co + Ceos (coti-*— <p) 4- a (30) 

where (i == 1, . . . , A:) are the time instants 

at which the measurements Zi were taken. Co 
denotes the mean level of Zi, C the amplitude, 
o) 2 = 27t/T the known angular frequency, 
and <p the so-called acrophase, that is, the 
phase at which the peak of the quantity z 


am indebted to F. Halberg, University of 
Minnesota, Minneapolis, for the peniussion'^ to use 
his data. 
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FIGURE 17. Application of the Cosinor method to telemetric measurements of intraperitoneal 
temperatures of adult female rat. Rat is subjected to special lighting regimen which influence 
the acrophase. Data were obtained in preparation for a space shot. 


occurs theoretically. Finally is the error 
term. 

In order to apply the standard technique 
of least squares estimation we assume that 
the errors £i are independently distributed 
normal variates with mean zero and common 
variance o^. 


The model is not linear in the unknown 
parameter but it can be linearized by 
rewriting 

C cos ( (oti — <p) — C cos 0)ti cos <p 

+ C sin o)ti sin (p (31) 
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and by substituting 

C cos <p^=x, C sin <p — y (32) 

Thus G and <p are replaced by new pa- 
rameters X and y. Notice that x and y can be 
interpreted as rectangular components of a 
vector of length C and polar angle (p. 

It is not the place here to present the 
method of least squares estimation. We men- 
tion only tha^we get minimum variance esti- 
mates X, % Co for the unknown parameters 
and an estimate for the unknown variance 
cr^ of the error term. 

Estimates for C and cp are found by solv- 
ing the equations (32) with x and y replaced 
by X and y^ respectively. 

If measurements are not only taken from 
one individual but from n comparable indi- 
viduals independently, we obtain a sample 
of estimates 

(33) 

The biological interest usually focuses on 
the amplitude and the acrophase. It is then 
natural to plot the n vectors with polar co- 
ordinates Cj‘, (pj or corresponding rectangular 
coordinates Xj, yj in an A:y-coordinate system. 
If our bivariate sample is roughly unimodal, 
it makes sense to calculate the components of 
the sample mean vector: 

x = (34) 

This mean vector in turn defines a mean 
amplitude and a mean acrophase for the 
sample. 

Let /Xaj and fjLy denote the rectangular 
coordinates of the true but unknown mean 
vector. We may then be interested in a con- 
fidence ellipse that covers the point {/Xa., fJLy ) . 
For this purpose we follow again a standard 


technique and determine the sample corre- 
lation coefficient r and the standard errors 
Sx and Sy of ^ and Then the point {fix, 
juLy) satisfies the inequality 

(35) 

Here F 2 , n— 2 denotes Fisher’s F for a preas- 
signed confidence coefficient. The inequality 
(35) defines the desired confidence ellipse. 

The confidence ellipse informs the re- 
search worker to what extent he can rely on 
the estimated mean amplitude and mean 
acrophase. If, for instance, the confidence 
ellipse does not cover the point (0,0), we 
conclude statistically that the amplitude is 
significantly different from zero. 

Conversely if the confidence ellipse covers 
the point {0, 0), then we have no reason to 
assume that there exists a nonzero amplitude. 
In this case, periodicity of the quantity z can- 
not be established. 

This method and related techniques, pro- 
posed by Halberg and coworkers, have be- 
come known as Gosinor method. For a more 
detailed account see reference 27. It is ex- 
pected that the Gosinor method will not only 
be useful in the area of biological rhythms 
as hitherto but also in some problems of 
orientation. 

There are other statistical methods which 
also deal with biological rhythms. For a re- 
markably clear treatment of periodic regres- 
sion, see reference 28 (chapter 17). 

DISCUSSION 

Williams : A number of us have been going to 
extremes to get not only the point on the circle or 
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vanishing point but also to develop a statistic that 
will express the degree of orientation of a track. 
The problem essentially is this : Given two sets of 
radiating lines, radiating in some primary direction, 
other than just plotting where they intersect a cir- 
cle, how can one determine whether these represent 
two different populations? Secondly, how can we 
develop a statistic that will express the straightness 
of the plot? 

BATsciiELEt : To express the degVee of orienta- 
tion we may subdivide the track into a moderate 
number of sections, each of them with the same 
flying time. We repIacO each section of the track by 
a straight line and Represent their directions by unit 
vectors. Then we calculate the mean vector as 
usual. Its angl€ serves as mean angle, and its length 
r is a measure of concentration. We may also use 
the angular deviation s [2 (1— r)]^ as a meas- 

ure of dispersion. 

Since each radiating line can be replaced by a 
straight line with a certain mean angle, a set of ra- 
diating lines can be considered as a circular sample. 
If such a set is a random sample from a population, 
circular tests can be applied. Likewise for two sets, 
circular two-sample tests are applicable. 

The straightness of a plot can be statistically ex- 
pressed in many different ways. One such way is 
using the length of the mean vector as described 
before. Another statistic or index would be the 
maximum deviation from a straight line joining ini- 
tial and end points of a track. Whatever statistic 
we use, we have to keep in mind that we lose some 
information and that there will be no definition 
that serves all purposes. 

Special difficulties arise if the animal follows a 
geographic line or is watching for a food source. 
Under variable atmospheric conditions the track 
could also be the result of a learning process. 

, Carr: What dp you think about that series of 
lO-miriute duty-cycle heading reports mentioned by 
Mr, Baldwin this morning? You suggested that 
these could not be evaluated statistically because 
the animal had perhaps learhed something between 
each successive heading. I don^t understand this. 
Why can’t each reported heading be considered in 
the same light as a single sally of an experimental 
bird in an orientation arena and regular circular 
statistics applied ? 

Batschelet: A statistic for the “straightness” 
of a track could be c|uite suitable to express how 
well an animal is keeping a certain direction. How- 
ever, the question may be more demanding : Is the 
animal heading for a particular destination or can 


the track be explained otherwise, say by random 
movements? With one single animal this, question 
cannot be solved by today’s statistical methods. We 
should first create appropriate models for animal 
migration. Such models would belong to the large 
area of stochastic processes. Then statistical meth- 
ods would have to he developed to test whether an 
observed track fits the model. A random sample of 
a few animals seems to be the easiest method to 
solve the problem of how well the animals are mo- 
tivated. 

Waterman: We are all very Concerned with 
these problems of experimental design and data 
analysis for orientation. One great difficulty is that 
we know almost nothing about the dyn^ics of ori- 
entation, consequently we do not really know what 
endpoints to use. We can take as many observations 
ad lib, but are they really independent? 

Another problem which is prominent in our re- 
search on polarotaxis is multiple peak orientation. 
It is difficult to decide how many significant orien- 
tation peaks you really have. In addition, we are 
also interested in knowing the location of these 
peaks and estimating the significance of their dif- 
ferences. , 

Batschelet : Multiple peak orientation occurs 
quite frequently. It is essential in each experiment 
to find out the reasons for this behavior before sta- 
tistical methods are applied. Once the theoretical 
peak directions are defined, statistics is able to test 
the goodness of fit or to decompose a multimodal 
distribution into unimodal distributions. However, 
counting peaks without a biological or physical 
model is a hopeless statistical enterprise. 

Carr: Suppose you put out a mechanical turtle 
with a built-in orientation system and yoii want to 
answer the question, Is the machine orienting or 
travelling at random ? Is there no statistical way for 
evaluating its capacity simply to hold a course? Is 
there no test for significance for a ^gment of a 
course, in terms of its adherence to a straightline or 
to some regular modification of a straight line? 

Batschelet: In principle it can be solved, but 
in the graphs I have seen, it cannot be handled 
that way. 

Baldwin : What is your opinion of the validity 
of comparing the paths of two turtles released si- 
multaneously, assumihg that thisre is no communi- . 
cation between them? . ’ 

Batschelet: Whether two turtles released at 
the same time move independently of each other 
cannot be tested by statistical methods. 
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